Strong and Electromagnetic Decays of X(1835) as a Baryonium State 
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With the assumption that the recently observed X(1835) is a baryonium state we have 
studied the strong decays of X(1835) — + yy'^'^Tr+Tr^, r^'^'^Tr'^Tr" and the electromagnetic decay 
of X(1835) — !■ 27 in the framework of effective Lagrangian formalism. In the present in- 
vestigation we have included the contributions from the iso-singlet light scalar resonances 
but we have not included the isospin violating effect. Our result for the strong decay of 
X(1835) ryV+TT^ is smaller than the observed data. The decay width for the radiative 
decay of X(1835) ~> 27 is consistent with the assumption that it decays through the glue- 
ball. In addition, the width for the strong decay of X(1835) TyTr+vr^ is larger than that of 
the strong decay of X(1835) rfn^TT^ due to the large phase space and coupling constant 
dNNri- From our investigation, it is not possible to interpret X(1835) as a baryonium. 

PACS numbers: 13.25.Jx,12.39.Mk, 13.40.Hq 

I. INTRODUCTION 

In 2005, the BES collaboration announced the observation of a resonant state termed X(1835) 
in the reaction J/tp ^ jX, X tI'k^'k^ [1]. A fit to this resonance with the Breit-Wigner function 
yields the quantum number J^^ = ^ and mass Mx = (1833.7 it 6.1 it 2.7) MeV, width Tx = 
(67.7 ± 20.3 ± 7.7) MeV and the product branching fraction Br(J/'0 7X(1835))Br(X(1835) 
r/'vr+vr") = (2.2ib0.4(stat.)ib0.4(syst.)) x 10"^. Actually, without include the final state interaction, 
the parameters of this resonance have been fitted to be Mx — 1859l'^o(stat)l'l25(syst) MeV and the 
total width F < 30 MeV in Ref . 4] . Since the discovery of the X(1835) state, many models have 
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been proposed to explain its properties 
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In the previous works, the X(1835) state has been conjectured to be a baryonium state 



8, 



18 
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pseudoscalar glueball state 
19[. Although there are many speculations, none of the above claims can be either 
confirmed or ruled out by the present experiments. In our present work, we have calculated the 
strong decays of X(1835) — > r]^'^TT^TT^ jIJ^'^tt'^tt^ and radiative decay of X(1835) — > 27 using the 
effective Lagrangian formalism by treating the X(1835) as a NN baryonium. This seems to be a 
reasonable approximation if one only considers the fact that the mass of the X(1835) is bit lower 



than the threshold energy of pp and nn (about 40 MeV). Our philosophy is that, assuming the 
X(1835) as a baryonium, if we can get the numerical results agree with the observed data the 
baryonium assumption is reasonable otherwise the baryonium picture can be ruled out, at least in 
this framework. The coupling of the X(1835) to its constituents can be described by the effective 
Lagrangian. The corresponding effective coupling constant is determined by the compositeness 
condition Z = which was earlier used by nuclear physicists 



used by particle physicists (see the references in 



study the newly observed charmed mesons 
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22( 1 and is being widely 



We had applied the above method to 



25l | and their decay properties which we had 



obtained agreed with the observed data. We had also employed the above technique to predict 

n 

the decay properties of the bottom-strange mesons [26]. In our present work, we have used a 
typical scale parameter Ax to describe the finite size of the baryonium. The value of Ax is fixed 
by considering the coupling constant is expected to be stable. For other interactions, we have 
used the phenomenological Lagrangian and have borrowed the relevant coupling constants from 
the existing literature. Using the above phenomenological approaches, we have analyzed the strong 
decays of X(1835) — > rj^'^ tt~^ 7r~ ,rj^'^Tr^Tr^ and radiative decay of X(1835) — > 2j. The result of the 
decay width of A'(1835) r]'Tr^iT^ is much smaller than the observed data hence the X(1835) 
cannot be treated baryonium. 

The paper is organized in the following way: In Section |IT1 we have calculated the effective 
coupling constant g^ using the compositeness condition and have discussed the effective Lagrangian 
formalism employed in our calculation. In Section IIIII we have calculated the strong decay widths 
of A'(1835) 7/'-'^7r+7r~ and radiative decay width of X(1835) 27 using the effective coupling 
constant g^^ and effective Lagrangian proposed in Section [III In section IIVI the important results 
and conclusions have been given. 



II. THEORETICAL FRAMEWORK 
A. Baryonium structure of the X(1835) state 

In this section we give the formulation for the study of the X(1835) as a baryonium state which 
can be thought of as a pp{nn) bound state. As stated earlier, the mass of the X(1835) is around 
40 MeV less than the threshold of pp{nn). The quantum number of the X(1835) is assigned to 
be J^^ = 0"+, and its mass is predicted to be nix = 1833.7 MeV [l|]. The effective Lagrangian 
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describing the interaction between the X(1835) and its constituents is given by 

^x{iS35)ix) = igx^i^) j dy^xiy^)Nix + ^y)-/5N{x-^y), (1) 
where the baryon doublet N is defined as 

N = 

The correlation function characterizes the finite size of the X(1835) as a NN bound state and 
depends on the relative Jacobi coordinates y and x. In the limit — ^ <5^(y)) the interaction 
given by Eq. ([1]) becomes local. The Fourier transform of the correlation function is 




^x{y 



(2vr)4 

In following calculation, an explicit form of has been used. The choice of ^x should be such 
that it falls off sufficiently fast in the ultraviolet region of Euclidean space to render the Feynman 
diagrams finite in the UV region. In this sense, one can also regard ^x as a regulator for the loop 
integral. In our work, we have chosen the Gaussian form for ^x 

^x{pI) = exp(-p|/A^) 

where j?_e is the Euclidean Jacobi momentum. Here Ax is a size parameter which parameterizes 
the distribution of N and N baryons inside the X(1835) baryonium. 



The coupling constant is determined by the compositeness condition 
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221 which 



implies that the renormalization constant of the hadron wave function is set to zero 



Zx 



T,'v{m 



2 ^ 

XV'Xj 



0. 



(2) 



Here, S^(m^) = (7^H^(m^) is the derivative of the mass operator which is represented by 
the diagrams in Fig. [1] given below. 





FIG. 1: Mass operator of X(1835) as a baryonium state. 
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The compositeness condition can be understood in the following: The renormalization constant 

1 /2 

Z-^ can be interpreted as the matrix element between the physical state X(1835) and correspond- 
ing bare state Xo(1835), i.e., (0|Xo(1835)|X(1835)) = z]/^(0|X(1835)|X(1835)) = so that 
Zx = means that the physical state should not be a function of the corresponding bare state 
which means that the physical state is a bound state. In our present work, the X(1835) is a bound 
state of pp{nn). In this sense, the compositeness condition excludes the possibility of the processes 
involving the X(1835) as an initial or a final state since each external X(1835) contributes a factor 

1/2 

Z^ to the relevant matrix elements. In addition, because of the interaction between the X(1835) 
and its constituents, the mass and wave function of the X(1835) have to be renormalized. 
Following Eq. ([2]) the coupling constant can be expressed as 



, „ , da dx 
./o 7o 



(1 + a)2 dfi'j^ 



(3) 



where a and x are both Feynman parameters and 
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am. 
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4(1 + q) 
a + ia'^xil - 
4(1 + a) 



-mx 



ml /A 



X 



X and a are both Feynman parameters. In deriving the expression ([3]), we have ignored the 
mass difference between proton and neutron and expressed the coupling constant in terms 
of the proton mass. To get the numerical result of g^^, we use mx = 1833.7 MeV jl|, mp = 
938.272 MeV [27] and vary the scale parameter from 1.0 GeV to 5.0 GeV. In Fig. [2] we show 
the Ax dependence of the effective coupling constant g^ . 




FIG. 2: The Ax dependence of the coupling constant g 
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Concerning that g^^ is expected to be stable against Ax, we choose the region of Ax as 
2.0 GeV < Ax < 3.0 GeV and get the couphng constant to be in the range = 2.37 — 2.55. 
Comparing our present result with that given in Ref. [3l where this coupling constant was esti- 
mated from experimental branching ratio of the X(1835) to pp decay in radiative decay of J/V' 
(by considering that X ^ pp occurs from the tail of its mass distribution and the value was found 
to be, \gxpp\ — 3.5), we conclude our result agrees with the result given there. In fact, using 
BR{X pp) ~ (0.04 - 0.14) assuming Tx < 30Mev\^ that Ref. adopte(Lone can get 

9xpp ~ 2-2 — 4.1. In addition, our conclusion is also consistent with that of Ref. [l5| which was 
based on the glueball assumption. Expressing the coupling constant gxpp in terms of g^^ which is 
the coupling constant between the X(1835) and glueball, one can get g^pp = 2.47 — 4.67. 

B. Effective Lagrangian for strong and electromagnetic decays of X(1835) 

In this section, we have discussed the effective Lagrangian for the calculation of the strong decays 
of X(1835) 7/('V+7r^ and electromagnetic decay of X(1835) — > 2j. The effective lagrangian can 
be divided into two parts, the free part £frcc and the interaction part £int- It should be noted 
that the electromagnetic interaction can be obtained by the minimal substitution (i.e., replacing 
the derivative operator of the charged particle with the covariant one = — ieQA^ with 
Q as the charge of the relevant particle). For the free Lagrangian, it involves states with quantum 
numbers = ^^,0^,0^ and 1^. 

■^free = -^free ~^ ■^free ~^ ■^free ~^ ^iree 

where 

= N{i^-mN)N 

J^Le = -^^(1835)(n + m2,)X(1835) 

-^7f(x)(n + ml)TT{x) - ^??(x)(n + m^,)r/(x) - ^r?'(x)(n + ml,)7]'{x) 

J^Le = -\a{D + ml)a-^fo{n + ml)fo 
1 



with Ff^i, = d^Au — dyA^ as the field tensor of photon and □ = d^d^. For computing the decays 
of the X(1835), we have treated the masses of proton and neutron and the masses of the triplet 
pions to be the same [1|, |27|] 

mx = 1833.7 MeV; = m„ = 938.27203 MeV 



m^o = = 139.57018 MeV; = 547.51 MeV; m^/ = 957.78 MeV; (4) 
while for the masses and widths of scalar mesons, we have adopted 28 1 



rua = 550 MeV; = 370 MeV; ruf^ = 980 MeV; Tf^= 64.6 MeV 

In the following calculation, we have included the finite width effects of the scalar mesons, that 
is, we have written the scalar meson propagators as 

Ds{k) 



k"^ — m| + imsTs 

The interaction Lagrangian £int used in our calculation has two parts, the strong part and 
the electromagnetic part 

n rstr , rem 

''-int — ''-int "T ''-int 

For the strong interaction Lagrangian we have C^ji^x (X-nucleon-nucleon interaction), C^fjp 
(pseudoscalar-nucleon-nucleon interaction), C^fqg (scalar- nucleon-nucleon interaction) and Cspp 
(scalar-pseudoscalar-pseudoscalar interaction) 

■^int = ^NNX + ^NNP + ^NNS + ^SPP 

The effective Lagrangian Cj^f^x given in Eq. ([T|) and C^^f^p and Cj^f^g can be expressed as 

^NNP = ^5AfiV7r%M75TA^5^7r + ^9^^^iV7^75Afa^r/ + ^5^^^,iV7^75A^a;,r/' (5) 
^NNS = 9nns^^S (6) 
J^spp = ■ ^/^^ ~ -^^fodi,^ ■ di^^ (7) 

where S is the scalar meson (cj and /o in our problem) and vr is the pseudoscalar meson matrix 



vr 



vr = ^ttV* = 

i=l \ V^TI" — 



V27T- — 



The coupling constants ff^Tv-n-' 9NNri 9NNr]' were determined via the J ftp hadronic decay 

while Pmns yielded by fitting the theoretical results of NN scattering with the observed 



data 



3ii 



(<7^jvJ7(4vr) ^ 14.8 
i9NN^/9NNn f ^ 3.90625 x lO'^ 

i9NNr,'/9NNn? - 2.5 X IQ-^ 



{dNNsfli^^) ^ 5.69 



The scalar-pseudoscalar-coupling constant 75pp was given in Ref. [28|] 



7^^^ = 7.27 GeV-i; 7/^^^ = 1.47 GeV"! 

For the electromagnetic interaction Lagrangian used in our calculation, it has two parts: 
(i) is from the gauge of the charged free nucleon Lagrangian, and (ii) is from the gauge of the 
nonlocal interaction 

rem _ f^emli) ^em(ii) 



where 



^Z^^ = eA,N^,^^N (8) 



where the Wilson line I{x,y,P) is defined as 

Iix,y;P) = r dz^A^'iz) 
To derive the Feynman rules for photons, we require the derivative of I{x,y;P). For this we 



33( 1 which implies that the 



have used the path-independent prescription as suggested in Ref. |32l . 
derivative of I{x, y; P) does not depend on the path P originally used in the definition. Also in our 
calculation of X(1835) 27, in principle we should expand the above expression to the second 
order but the diagram with photons from this vertex does not contribute since the trace of gamma 
matrices vanish. 

III. STRONG AND ELECTROMAGNETIC DECAYS 

Having discussed the effective coupling constant and the effective Lagrangian, we are in the 
position to calculate the decay properties of the X(1835). In this section, we have calculated the 
strong decays of X(1835) — > ■r]^'\^T:~ and also the radiative decay of X(1835) — > 27. 

A. Strong decays of X(1835) Ty^'^Tr+Tr" 

For the strong decays of X(1835) Ty^'^vr+vr^, the Feynman diagrams of Fig. [3] and Fig. U] 
contribute. All the diagrams listed in Fig. [3] are from the one-pseudoscalar meson- nucleon- nucleon 
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vertex while the diagrams hsted in Fig. H] are from the scalar resonance contributions. For the 
isospin symmetric case following relations among matrix elements exist 




(£>) {E) (F) 

FIG. 3: Diagrams contributing to the strong decay of decay of X(1835) — > jf-'^-K^-K^ without scalar resonance 
contribution. 




FIG. 4: Diagrams contributing to the strong decay of decay of X(1835) jf-'^-K^-K witli scalar resonance 
contribution. 

It should be noted that to include the isospin violating effect, the diagrams in Fig. [5] and Fig. [6] 
should also be considered. For isospin symmetric case the matrix elements for the diagrams of 
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Fig. [5] and Fig. [6] have the following relations 

In our present work we have considered isospin symmetric case and hence diagrams of Fig. [5] and 
Fig. [6] do not contribute. In addition, the diagrams with p^r/^'^d^TT^ vertex also have not been 
considered due to the G-parity conservation. 



7,0 




FIG. 5: Diagrams contributing to the strong decay X(1835) — ^ rj'^'^'K^Tr from the p meson exchange. 




FIG. 6: Diagrams contributing to the strong decay X(1835) — * Ty'^'V+Tr from the two- meson- nucleon vertex. 

In the following calculation, we label the momenta of the relevant particles according to the 
scheme X{p) —>■ vr+(gi) -|- vr~(g2) + ii^'\q3)- The partial decay width is related to the invariant 
matrix element M(p ^ qi + q2 + Q3) by the relation 

r(X(1835) ^ ^+7r-r/W) = [ \M^\d<^ 

2mx J 

9 



where d$ is the Lorentz invariant phase space volume element 

3 3 



3 



with E\ = \J m\ + , E2 = \J m\ + and £^3 = ^rv?^^^,^ + q| . After integrating the delta 
function over the solid-angle elements d£l\ and d£l2 and treating the X(1835) as an unpolarized 
particle, the partial decay width can be expressed as a two dimensional integral 

r(X(1835) ^ r/('V+^-) = —4 / \M'^\dEidE2 (10) 

647r'^mx J 

The matrix elements are calculated by evaluating the loop integral. For example, the matrix 
element M^^-j for the corresponding diagram (A) in Fig. [3] is 

TV{(|^ + m)75[(^ -i>)+ m]^375[(^ - ^1 - ^2) + m]^275[(^ - ^1) + rn]^i75} 
(/c^ — m^) [(A; — p)^ — m'^] [{k — qi — q2)'^ — m"^] [{k — qi)"^ — inn?\ 

where Qeff = Oxd^qfj-j^dNNr]'^')- After performing the trace calculation, the matrix element can be 
decomposed in terms of the tensor structure 

where a's are functions of the external momenta and D^s are the loop integrals. Their explicit 
forms are given in Appendix [Xl 

The results for the decay widths of T{X r/^'^vr+vr^) in the energy region A^ = 2.0 — 3.0 GeV 

are 

T{X r/'vr+Tr") = 0.580 - 1.273 MeV; T^{X Vvr+vr") = 0.335 - 0.400 MeV 
T{X ?77r+7r") = 6.522 - 13.29 MeV; T^{X r/7r+7r~) = 1.550 - 1.926 MeV 

where the upper index P means that the results are from the pure pseudoscalar processes illustrated 
in Fig. [3l The above decay widths increase with increase in Ax. To obtain the above results, the 
coupling constant calculated before and the coupling constants given above were used. Using 
the central value of the total width r(X(1835)) = 67.7 MeV [l{], the branching ratios turn out to 
be 

BR{X r/'vr+TT-) ~ 8.57 x 10"^ - 1.88 x 10"^ 
BR(A: ^ ryvr+vr") ~ 9.63 x 10^^ - 1.96 x 10"^ 

10 



Using the result BR(J/?/' ^ jX) ~ (0.5 - 2) x lO'^ 
fraction is obtained 



, the following product for branching 



BR{J/ip jX)BR{X -q'-K+Tr-) ~ (0.428 - 1.714) x 10"^ - (0.94 - 3.76) x lO^'^ 

which is much smaher than the observed data. The uncertainties in the parentheses are from the 
uncertainty of BR( J/-^ iX). The large uncertainty comes from the measurement of BR(J/'0 
'jX). In addition, the product of branching fraction BR(J/?/; jX)JiIi{X rjir^TT^) yields 

BR(J/V' 7X)BR(X r/TT+vr^) ~ (0.418 - 1.926) x 10"^ - (0.963 - 3.852) x lO"'^ 

where the uncertainties in the parentheses are also from the uncertainty of BR(J/'0 — > ^X). 

Our calculation shows that the strong decay width r(X(1835) — > r/'vr'^vr^) based on the bary- 
onium assumption in the energy scale 2.0 GeV < Ax < 3.0 GeV is much smaller than the 
data which leads to the conclusion that the X(1835) may not be a baryonium. In addition, 
we have also predicted the strong decay width of T{X ryvr^vr^) should be larger than that 
of r(X(1835) r/'vr+vr^) if the X(1835) is a baryonium due to the large phase space and coupling 

constant ^Ar^Vr;- 

B. Radiative decay of X(1835) 2-f. 

The X(1835) state can decay into two photons. Since the X(1835) state is a pseudoscalar state 
the radiative decay is an anomalous process. The matrix element can be written as 



iM'^'iX (1835) ^ 27) = aeme^ua/3Paqi3Gx^ 



77 



where q and p are the momenta of the two final photons. Using the above expression the decay 
width is given by 

r(X(1835) ^ 27) = ^iMpM^I^ = -^aLm^G^,, 
Qirnix rnx oZn ' ' 

where |pcm| = ^T^x/2 is the three-momentum of the decay products. 

In our present model, the decay X(1835) — > 27 happens via the process given by the diagrams 

in Fig. [71 Diagrams (A), (B) and their corresponding cross diagrams arise from the gauge of the 

nonlocal interaction ([9]). Diagram {A) and its cross one are from quadratic terms of A^ in the 

expansion of Eq. ^ while diagram (B) and its cross one are from the linear terms of and the 

gauge of the proton free Lagrangian ([8]). Diagrams (C) and (D) arise from Lagrangian given by 
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FIG. 7: Diagrams contributing to the radiative decay X(1835) 2j. 

From our analysis neither diagram (^4) nor diagram (B) contributes to the total matrix element 
due to the vanishing of the trace of gamma matrices. So we need to calculate only the diagrams (C) 
and (D) which are the same as that calculated in the triangle anomaly problem. Since the discovery 



of the triangle anomaly 



35 



361 ] ■ the calculation of these diagrams have been discussed widely in the 

literature. We had discussed the ambiguities in the calculations induced by regularization, Dirac 

trace, and momentum shifts [3^]. From our calculation 

2m f°° 1 

<^X77 = fx— rf / daida2da3—— exp{arg^J Aj^} 

" 71" A Jo (1 + ai + 02 + as)^ 

where the upper index NL corresponds to the Nonlocal case and 

Ill 1 ^ 

argem = (tt + a2)(- + a2 + a3)"ix + (t + "2)"^^ - Vttim^ 

1 + Li=iai ^2 4 .^^ 

Using the values of the parameters we present the numerical results now. For the effective 

coupling G^^, and for the scale parameter in the range 2.0 GeV < Ax < 3.0 GeV, we get the 

result 

G^^^^ = 0.6813 GeV"^ - 0.3804 GeV"^ 
and the corresponding electromagnetic decay width 

r^^;^ = 1.516 KeV - 0.4726 KeV 
12 



9x 


^ X — >r/'7r+7r~ 








^X^2i 


2.55 - 2.37 


0.580 - 1.273 


6.522 - 13.29 


0.290 - 0.637 


3.261 - 6.645 


1.516 - 0.4726 



TABLE I: The Ax dependence of effective coupling constant and decay widths in the region 2.0 GeV < 
Ax < 3.0 GeV(The strong decay width is expressed in unit of MeV while the electromagnetic decay width 
is expressed in unit of KeV). 



Both Gx77 '^^■y decrease with increase in Ax- 



X77 



The radiative decay has been investigated in Ref. 



3 



treating the X(1835) glueball. The 



result obtained for the decay width Tx'yy = 1.1(0.31 — 1.1) KeV agrees with our result. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this work, the strong decays of X(1835) rj^'^ir^Tr^ and electromagnetic decay of X(1835) 
27 have been calculated using the effective Lagrangian method. In our work we have treated the 
X(1835) as a baryonium. To fix the only free parameter A^ we postulated that the coupling 
constant has to be stable against Ax- With this assumption, we varied Ax from 2.0 GeV 
to 3.0 GeV. In the above region the strong decay width of X(1835) — > rj'ir'^Tr' is much smaller 
than the observed data but our prediction of the electromagnetic decay width of X(1835) 27 
is in agreement with the result where X(1835) decays through glueball. In addition, we have also 
calculated the strong decay width r(X(1835) — > r/Tr^vr") explicitly. The calculated width is much 
larger than the partial width of r(X(1835) r/'Tr+vr") which is consistent with the direct analysis 
of the phase space and the coupling constant. 

In the baryonium picture, other decay modes of X(1835) can also be calculated. Using the 
isospin relation we get 

r(X(1835) ^ ?7'7r°7r°) = ^T{X{1835) ti'tt+tt-) = 0.290 - 0.637 MeV; 
r(X(1835) r/7r°7r°) = ^r(X(1835) 7?7r+7r") = 3.261 - 6.645 MeV 

The other three-pseudoscalar strong decay channels are either isospin symmetry violating processes 
(tt+tt^tt'^ and 37r'^) or OZI rule suppressed (with Kaon meson in the final state). The four strong 
decay channels discussed above are dominant among all the three-pseudoscalar channels. We have 
listed the effective coupling constant and their decay widths in the region 2.0GeV < Ax < 
3.0 GeV in Table. H 

It should be noticed that in principal, the finite width effect should be included by introducing 
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the Breit-Wigner distribution function. However, this wih suppress our results and our final con- 
clusion will not be changed. Moreover, there are also uncertainties from the sigma meson mass and 
width. Here, we applied the results yielded by unitarizing the vrvr and ttK scattering amplitudes. 

To conclude, we have studied the three-pseudoscalar meson and two-photon decays of X(1835). 
The strong decay width r(X(1835) rj'ir'^Tr^) is smaller than the experimental data while the 
two-photon decay width agrees with the result where X(1835) was assumed to decay via the glue- 
ball assumption. From our results X(1835) cannot be treated as a baryonium, at least in the 
framework of the composite model as applied in this paper. We have obtained other dominant 
three-pseudoscalar meson decay channels from the isospin relations. To confirm the structure of 
X(1835) further theoretical analysis is necessary. 

APPENDIX A: DECOMPOSITION OF ONE LOOP INTEGRAL. 



For the one loop integral of diagram (A) of Fig. [3l after performing the trace calculation we get 
the following decomposition 

2 r ^4/ 



iM, 



(A) 



where 



a 



a 



a 



. TV{(^ + m)75 [(^ -^)+ m]^375[(^ ■ 



(k"^ — m?) [{k — p)"^ — m^] [{k — qi — q2)'^ — m"^] [{k — qi)'^ — nn?\ 



(2m)2 



m 



+m 



m 



+m 



m 



+m 



2m 



P ■ qm ■q2-p- q2qi ■ qs - 2gi • 92^1 -qs-qi- + p • Q1Q2 • qs + 2gf ^2 • qs 



■ p ■ q-m ■ q2ql - p ■ 92^1 • q?,ql - p ■ qmlql + p ■ Q1Q2 ■ qml 



qi ■ ^2^3 - 2^2 • qsqi 
- 2p • q2qlq3 + 4p • qsqjq!^ + 4p • qiqi ■ 9293 - 2p • 93^1 • 9291 + 2p • ^291 • gs^i 
-4qi ■ 5291 • qsP^ + 2p • 919^93 - 2qi ■ q^qlp^ - 2p ■ qiq2 ■ q^i + 2q2 • qsqfp^ 
4p ■ qiqt^q^ + 4p ■ 529193 " 4p • 939^92 + ^Qi ' I^P^ll - ^92 • 93P''9i 



3p • 9391 • 92 - 3p • 5291 • 93 + 291 • 9291 • 93 + 9i • 9392 + 3p • 9i92 • 93 - 29^92 • 93 



92 • 939i - 91 • 9293 
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and 



(27r)4 2''(A;2-m2)[(/c-p)2-m2]p-gi-g2)2-m2][(A;-gi)2-m2] 
It is to be noted that due to the relation 

k-p = -H[{k-pf -m'^]-{k'^ -m'^)-p^ 



A[{k-qi? -m']-{k' -m')-ql 



k ■ qi 



k-q3 = -^^[ik-pf-m'^]-[{k-qi-q2f-m'^] + {qi + q2f 

the above vector, two- and three- rank four-point integrals can be expressed in terms of scalar 
four-point and three-point integrals 



+m 



Po^Dq + /3i23C'i23;0 + /?124C'l24;0 + /5234C'234;0 



dT2 



T2 . 



oT2 , 



92 • 93^1 - qi ■ 9253 



q2 ■ 93[C'l23;0 " C'234;o] " 9l ' 92[C'l34;0 " C'l24;o] 

-{q2 ■ qml + qi ■ 52 [(gi + q2f - P^]}Do 



with 



+m 



qi ■ q2[{qi + q2f - q'^] + 2^2 • 93^1 

- 2g • q2ql[{qi + ^2)^ - q^] + 4g • qsqliqj - {qi + ^2)^] 

+Aq ■ qiqi ■ 52 [(gi + ^2)^ - q^] + 2g • gs^i • ^2?? - 2g • 9291 • 939i 

+4gi ■ q2qi ■ qzq^ + 2g ■ gigi[(gi -|- ^2)^ - q^] + 2gi • qmlq^ + 2? ■ q'iq'2 ■ qzql 

-2^2 • 93919^ 



/?234 = 2m'' 92 • 93 + 



/^134 



--<; m'^gi ■ q2 + m 



2g • 9391 • 92 - 2g ■ 9291 • 93 + 49i ■ 9291 • 93 + 291 • 9392 + 29 • 9i92 • 93 



-292 • 939i 



29 • 929i + 49 • 91 91 • 92 - 491 • 9291 • 93 + 29 • 9i92 - 29i • 9392 



+292 • 939i 



/^124 



— mrqi ■ q2 + m 



29 • 929i + 49 • 939i - 49 • 91 91 • 92 - 29 • 91 9I 
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^ I - 2m'^q2 -qs + m - Aq ■ q^qf - 2q ■ q^qi ■q2 + 2q- q2qi ■ qs - 2q ■ qiq2 ■ qs 
= 3p ■ qsqi ■ q2 - 3p ■ q2qi ■ qs + 2qi ■ q2qi ■ qs + qi ■ q^ql + 3p • qiq2 ■ qs - 2qlq2 ■ qs 



Pi23 = -7;{ ~ 2m^q2 -qs + m 



dTI 



Pi24;fi = "^P ■ Q2qi - 2p ■ qiq^ 
Pl34;tM = 2p • 9192 - 2P • 9291 - 291 • 9392 + 292 • 939l 
l3234;n = 291 • 9392 " 2p • 9392 - 292 • 939l 



= -^|29 • 9i[(9i + 92)' - 9'][-(9i + 92)' + 9?] + 29 • 929?[(9i + 92)' - 9'] 
-29 • 939? [-(91 + 92)^ + 9?] + 291 • 939^ [-(91 + 92)^ + qf] + 292 • 939^? 
= -U - 29 • 91 [(91 + 92)' - 9'] - 29 • 92 [(91 + 92)' - 9'] 



2 



+29 • 939i - 291 • 939^ - 292 • 939^ 



Pi24 = -2 ^ 29 • 91 [(91 + 92)^ - 9 ] - 29 • 939r + 291 • 939' 
PI34 = 29 •92K91+ 92)^ -91 + 292- 939' 



and 



W23;{0;m} j (27r)4^^ 2' ' {k'' - m?)[{k - pf - m?][{k - qi - q2f 

^ f d'^^ P\\ 

Ci24;{0;m} = J J^^^^^ ' ^ 



2 {k^ — m^)[{k — py — ni?][{k — qiy — m^] 
^234;{0;/.} J ^2tt)* ^^ 2 ^ ' [{k - p)^ - m^][{k - qi - q2y - m^][{k - qi)^ - m^] 



d^k p.. {l;k^} 



m'^)[{k — 9i — 92)^ — m2][(/c — 91)^ — m^] 

Using the above, the matrix element iM(^) can be expressed in terms of the scalar, and vector 
C and D functions. For the scalar, vector C and D functions one can evaluate the momentum 
integral explicitly and yield the following results. 

i 1 f°° 1 1 

^{°^'^> = 16^ Jo ^«i'^«2^«3da4^^^:^{l;3^Pi,;,}exp{arg^/A3,} 

i 1 f°° 1 1 

Ci23;{0;/.} = 'Jq^J^ daida2da^ ^ {1; 3^7^-^0123;^} exp{argc^^3/Al^} 

i 1 1"^ 1 1 

Ci24;{0;m} = ~I^A^ /, daida2da3 ^ {1; ^^7^^0124;/^} exp{argg^^^/A|} 

i 1 1 1 

C234;{0;m} = "le^^Aj daida2da^ ^ {1; ^7^^-^0234;/.} exp{argg^3^/A^} 



16 



i 1 f°° 1 1 

Ci34;{0;m} = 'J^JT daida2da3j^^-^^{l; Y^^Pc,^,■,^.} exp{avgc,J Aj,} 



where 



n 

an = y~^«i 

i=l 

,1 

,1 



PD■,^l = i-^ + (X2)Pf, + a3qi-f, + a4{qi + q2)fj, 



^Ci23;m = (2 +«2)P/i + a3(gi+g2! 



1 



^Ci34;m = + "2(^1 +^2)^ + "391; 



arg^ = - ^ ^ P|) + (7 + 02)^^ + asg'i + a^iqi + 0^2)^ - a4"i^ 
1 + a4 4 

a'^SCm = -j^^^Cm + (^ + "2^ + a3(9i + 92)^-a3m2 

argci24 = -Y^^^Ci24 + + "2^ + asqj - a^ 
1 1 

argc234 = 0,3 -^^234 + (4 + "1)^^ + "2(91 + (12? + a^ql - asm^ 

argcm = - -I -Pci34 + 7^^^ + "2(91 + 92? + asqj - asm^ 
i + aa 4 
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